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Abstract 

Here we study a noise induced transition when the system is driven by a noise 
source taken as colored and non-Gaussian. We show -using both, a theoretical 
approximation and numerical simulations- that there is a shift of the transition 
as the noise departs from the Gaussian behavior. Also, we confirm the reentrance 
effect found for colored Gaussian noise and show the behavior of the transition 
line in the phase-like diagram as the noise departs from Gaussianity in the large 
correlation time limit. 
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1 INTRODUCTION 



During the last three decades a wealth of research results on fluctuations or 
noise have lead us to the recognition that in many situations noise can ac- 
tually play a constructive role that induces new ordering phenomena. Some 
examples are stochastic resonance in zero-dimensional [1] and extended sys- 
tems [2,3], noise induced transitions [4], noise induced phase transitions [5,6], 
noise induced transport [7,8,9], noise sustained patterns [10], etc. 

Most of the studies on the noise induced phenomena indicated above assume 
that the noise source has a Gaussian distribution (either white or colored). In 
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addition to the intrinsic interest in the study of non-Gaussian noises, there are 
some experimental evidences, particularly in sensory and biological systems 
[11], indicating that at least in some of these phenomena the noise sources 
could be non-Gaussian. The use of non-Gaussian noises in studies on noise 
induced phenomena is scarce mainly due to the mathematical difficulties. This 
is in contrast with the existence of some analytical tools when working with 
Gaussian (particularly white) noises. 

Here we present some results on one of those noise induced phenomena when 
driven by a noise source taken as colored and non-Gaussian. It corresponds to 
a noise induced transition like those discussed in [4]. The problem we discuss 
here corresponds the so called genetic model which follows from a particular 
biological modeling or from a chemical reaction [4,12]. In Ref. [13] the effect of 
a colored Gaussian (Ornstein-Uhlenbeck) noise on that model was analyzed 
in detail and a novel reentrance phenomenon in the phase diagram was found. 
In order to study the effect that the non-Gaussian nature of the noise can 
have in the transition, we have used here a particular distribution [14] whose 
departure from the Gaussian behavior is governed by a parameter q: while 
q = 1 corresponds to Gaussian noise, the value q > 1 produces a long tail 
distribution and q < 1 gives a cut-off distribution. As discussed for other 
noise induced phenomena in [15], we have found that the departure from the 
Gaussian distribution for the noise produces a strong effect. In this case such 
an effect corresponds to an anticipation of the transition when the distribution 
has long tails, and to a retardation when it is cut-off. This effect could have 
interesting consequences in the case of noise-induced-phase transitions [5,6]. 

In the next section we briefly discuss the form and properties of the non- 
Gaussian noise source. After that, we present the model we analyze here and 
a simple analytical approximation to treat the non-Gaussian noise that allows 
us to determine some of the general features of the phase diagram. We con- 
tinue discussing numerical simulations that support and extend the theoretical 
results. In the last section we draw some conclusions. 



2 NON-GAUSSIAN NOISE AND ITS PROPERTIES 

We start considering the following general form of a Langevin equation 

x = f(x)+g{x)r j (t). (1) 



As described in the Introduction, we assume that the noise term n(t) has a 
non-Gaussian distribution. More precisely, we consider that rj (t) is a Marko- 
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vian process generated as the solution of the following Langevin equation[14] 

(2) 



■n = —— v M + - — £(*), 



r drj 



being £(£) a standard Gaussian white noise of zero mean and correlation 
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Although we believe that our results are quite general, such a particular form 
for the noise rj(t) allows us to easily control the departure from the Gaussian 
behavior by changing a single parameter q. D and r are noise parameters re- 
lated to the noise intensity and the correlation time, as we now detail. The 
stationary properties of the noise 77, including the time correlation function, 
have been studied in [16] and here we summarize the main results. The sta- 
tionary probability distribution is given by 

-1 
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where Z q is the normalization factor. This distribution can be normalized only 
for q < 3. The first moment, (77) = 0, is always equal to zero, and the second 
moment, 



is finite only for q < 5/3. Furthermore, for q < 1, the distribution has a cut-off 
and it is only defined for \r/\ < r/ c = y / r ^ D <? - ) ■ Finally, the correlation time r q of 
the stationary regime of the process rj(t) diverges near q = 5/3 and it can be 
approximated over the whole range of values of q as r q m 2r/ (5 — 3g). Clearly, 
when q — > 1 we recover the limit of 77 being a Gaussian colored noise, namely 
the Ornstein-Uhlenbeck process, £, u{t), with correlations (C,ou{t)C,ou{t')) = 
R e -\t-t \/r anc [ probability distribution P s *(^ otl ) = Z~ l e~^D^° u . 

In Ref. [16], an effective Markovian approximation to the process 77 via a path 
integral procedure was obtained, that is analogous to the "unified colored noise 
approximation" obtained in [17] for the case of Ornstein-Uhlenbeck noise. 
Such an approximation allowed us to get quasi-analytical results for the mean- 
first-passage time or transition rate. These results and its dependence on the 
different parameters in the case of a double well potential were compared with 
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extensive numerical simulations with excellent agreement. That approximation 
was also exploited to study stochastic resonance with theoretical results [18] 
that are in very good agreement with experimental ones [19]. 



3 gsil APPROXIMATION FOR THE NON-GAUSSIAN NOISE 



The "effective Markovian approximation" introduced in Ref. [16] and indi- 
cated above has been shown to have some drawbacks. As discussed in [16], 
such an approximation renders a non normalizable stationary probability dis- 
tribution function, P^irj). This seems to be a general characteristic of such 
an approximation that, at least in the case discussed in [16], was solved intro- 
ducing an adequate and controlled cutoff. However, it is not clear that such 
a methodology could be always used, and particularly for the models we can 
analyze in relation to noise induced transitions. Hence, for the kind of models 
we use in order to study the effect of non-Gaussian noises on a noise induced 
transition we have resorted to a more simple, albeit physically reasonable, 
approximation. Our aim is at least to obtain some of the main trends in the 
region \q — 1| <C 1 (both for q < 1 and q > 1). The idea is that in such a region 
the non-Gaussian noise will only slightly depart from the Gaussian behavior. 
Hence in Eq. (2) we adopt 



1 dVgjrj) _ r] 
t drj t 
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= a(r)?7, 



5 — 3q 

with a(r) = — ; -. This implies a kind of renormalized Ornstein-Uhlenbeck 

2r(2 -q) 

(Gaussian) noise. Using this form of the noise we proceed to study the problem 
in Eq.(l). 

The genetic model is described by the following equation [4,12,13] 



x = — — x + \x{l — x) + x(l — x)r)(t), 



(7) 



and, according to the discussion above, we consider in this work that rj(t) is a 
non-Gaussian noise described by Eqs.(2-3). 

The analytical treatment simplifies if we make the change of variables 
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so that the original multiplicative stochastic differential equation (7) trans- 
forms into a new one but now with an additive noise 

y = -sinh(y) + A + rj(t). (9) 



From now on, as in Ref. [13], and in order to simplify the algebra, in what 
follows we adopt A = 0. 

The use of the approximation Eq. (6) (that we reiterate could be valid only 
for \q — 1| <C 1) together with the "unified colored noise approximation" [17], 
yields 

'"'WpiJ+^tt^ 1 ' 1 ' (10) 

L a(T) J L a(T) J 



with f(y) — — sinh(y). This result can be obtained in two ways: (a) Applying 
a direct adiabatic-like elimination procedure to the Langevin system given 
by Eq. (9) and Eqs. (2,6) [17]. That is, taking the derivative of Eq. (9) with 
respect to t and setting y = 0. (b) A formal one, through the application of the 
path-integral formalism to the indicated non Markovian Langevin equations, 
supplemented with a consistent Markovian approximation scheme [20]. 

Since, by the use of these "drastic" approximations, the problem has been re- 
duced to Eq. (10) which contains a white, Gaussian, noise, it is easy now, using 
standard techniques, to obtain the stationary probability density function. For 
the original variable x, it reads 

P st (x) = Nexp[F(x)}, (11) 



where iV is a normalization constant and 



F(x) =- 




:i2) 



+ ln 1 



-\n(x(l-x)). 



A noise induced transition [4] appears when this density function changes from 
being unimodal to bimodal. More specifically, the point x — 1/2 changes 
from being a maximum to being a minimum. The condition of extremum at 



5 



x — 1/2 and the condition for finding the boundary separating the regions 
where x = 1/2 is a maximum or a minimum yield 
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such that for a fixed r, the distribution is bimodal for D > D{t) and it is 
unimodal for D < D{r). In the limit q — > 1, this expression reduces to the 
result found in [13] for q = 1: 



Figure 1 plots the -D(r) boundary (the "transition line") for different values of 
q in the (D, r) plane. Several features can be noticed in this figure. It is clear 
that a minimum value of the noise intensity D is always needed to induce the 
transition. The dependence with the parameter q appears as an advancement 
of the transition to smaller values of D when q > 1 and a retardation when q < 
1. We believe that this will be a general behavior for long tail (power-law) and 
cut-off noise distributions, respectively. A closer look shows the persistence 
of the reentrance effect previously found in [13]: for some fixed values of the 
noise intensity D it is possible to obtain a transition, first from unimodal to 
bimodal and then back to unimodal, just by varying the time correlation r of 
the noise. 

These main features, as obtained from our very simplifying assumptions are 
found to be in qualitative agreement with the results coming from numerical 
simulations of the (not simplified) set of equations (2,3,7) (details of the nu- 
merical integration scheme are given in the appendix). In figure 2 we show 
that the probability distribution can switch from unimodal to bimodal just by 
increasing the parameter q, for fixed values of D and r. The reentrance effect 
is illustrated, for q — 1 in figure 3. 

The numerical simulations show that there is a crossing of the transition lines 
for q < 1 and q > 1, an aspect which is not reproduced by the approximated 
theoretical results. Another aspect observed in figure 1 is that for large r, the 
theoretically predicted transition line, has a smaller linear slope for increasing 
values of q. In fact, according to Eq. (13) we have D ~ 3(2-1) T * n li m h- 
There seems to be not such a strong q dependence of the slope of the curves 
for large r in the numerical results. Otherwise, the numerical and theoretical 
results are in good qualitative agreement. 
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Fig. 1. Phase diagram D vs. r for different values of q. Probability distributions 
are bimodal for points above the curves and unimodal below. Solid lines are the 
theoretical predictions, as given by Eq. (13) for q = 0.5, q = 1, q = 1.5 (from 
bottom to top). The numerical results are indicated by symbols joined by a dashed 
line for q = 0.5, a dotted line for q = 1, and a dash-dotted for q = 1.5. It is clear 
the reentrance effect and the shift of the transition line. 



We have presented some preliminary results on a noise induced phenomenon 
driven by a colored and non-Gaussian noise source that was generated by a 
g-distribution. The phenomenon studied here is a noise induced transition like 
those discussed in [4], and corresponds the so called genetic model. We have 
resorted to a simple approximation -that avoids the difficulties of the effective 
Markovian approximation in [16]- introduced with the idea that for \q— 1| <C 1, 
the non-Gaussian noise will only slightly depart from the Gaussian one, and 
its effect being only a correction to the Gaussian behavior. The resulting Eq. 
(13) shows some of the right tendencies with q > 1 or q < 1, although some 
features observed in the simulations (crossing of the transition lines and slight 
dependence of the linear slope for large r) are not reproduced by the simplified 
theoretical treatment. In addition, the stationary probability density function 
defined in Eqs.(ll) and (13) clearly shows, for constant values of D and r and 
as a function of q, the transition from unimodal to bimodal behavior. Again, 
these general results and tendencies, together with the reentrance effect, have 
been confirmed by extensive numerical simulations as shown in figures 1, 2 and 
3. However, the agreement between theory and simulations is only qualitative, 
as the first are unable to reproduce some of the details revealed by the second. 
For instance, from the numerical simulations a more complex behavior of the 



4 FINAL REMARKS 
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Fig. 2. Simulation results for the stationary probability distribution P s t{x) for fixed 
values D = 2, r = 0.4 and different values of q: q = 0.5 (*), q = 1 (x), q = 1.5 (+). 
It is clear that it is possible to induce a transition from unimodal to bimodal just by 
changing the parameter q, measuring the departure from the Gaussian distribution 
for the noise source. 

2.5 i . . . . . . 1 
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x 

Fig. 3. Simulation results for the stationary probability distribution P s t{x) for q = 1 
(Gaussian case) for fixed D = 1.95, and different values of r: r = 0. (x), r = 0.4 
(+), r = 1.0 (*). The reentrance effect indicated in the phase diagram D vs. r is 
apparent.. 
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transition line for large r is apparent, while it is not present in the theoretical 
predictions (see Fig. 1). 

These results are in the same line that a series of previous studies on noise 
induced phenomena where the noise source was non-Gaussian. In all the cases, 
we have found that the response of the system results to be enhanced or to 
have a strong dependence for values of the parameter q corresponding to a 
departure from Gaussian behavior. Some of such studies, including stochastic 
resonance [18,19], gated traps [21], Brownian motors [22], have been reviewed 
in [15]. 

Even though our results are so far only preliminary and further studies are 
required, we expect that the phenomenon discussed here will have remarkable 
effects on the noise induced phase transitions [6] as well as in the problem 
of coupled ratchets [8,9]. For instance, we can expect changes in the phase 
diagram of the first, or in the hysteresis cycle of the second, that could have 
interesting consequences in technological applications. Such studies will be the 
subject of further work. 
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APPENDIX: NUMERICAL SIMULATIONS 



In order to test our theoretical predictions, we have resorted to numerical sim- 
ulations of the system indicated by Eqs. (7) and (2,3). For the numerical 
integration of those equations we have used a second order stochastic Runge- 
Kutta type algorithm [23]. For the general problem 



x = f(x) + g(x)r](t) 

fl = q{v) + e£(t), (15) 

where £(t) is a white noise of mean equal to zero and correlations (£(£)£(£')) = 
5(t — t'), the algorithm is as following: time is discretized in steps of size h 



9 



and the value of the process rj(t) at the discretized times is obtained by the 
recursion relation 

k(t) = hq(r)(t)) 
l(t) = eh 1 ' 2 G(t) 

rj(t + h)= V (t) + £ [q(r,(t)) + q(r,(t) + k(t) + l(t))\ + l(t), (16) 

with the initial value 77(0) = 0. Here the numbers G(t) for different times t 
are independent Gaussian random numbers of mean equal to zero and vari- 
ance equal to one. They have been obtained by using a particularly efficient 
generator [24]. In the case q < 1, it might occur that the generated value of 
G(t) is such that r](t + h) is outside the cut-off interval (—r] c ,r] c ). If this hap- 
pens, that value is discarded and a new G(t) is generated until the condition 
\rj(t + h)\ < rj c is fulfilled. We have found that the percentage of discarded G(t) 
values is less than 5 x 1CT 4 in the worst case analyzed here corresponding to 
q = 0.5, t = 4. No such a problem ever occurs for q > 1. 

Once the values of rj(t) at the discretized times have been obtained, the inte- 
gration of (15) proceeds with a standard Runge-Kutta second order method: 

k x (t) = f(x(t)) 
x(t + h)=x(t) + ^[f(x(t)) + f(x(t) + k x (t))+ 

g(x(t))r)(t) + g{x{t) + k x {t)) V {t + h)}. (17) 

We have used a time step h = 0.01 and the histogram for P st (x) has been 
obtained from 10 8 values of x separated by a time t — 1 (corresponding, hence, 
to 100 integration steps). 
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